Less than 10 percent of star formation in z ~ 0.6 massive galaxies is 

triggered by major interactions 
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ABSTRACT 

Both observations and simulations show that major tidal interactions or mergers 
between gas-rich galaxies can lead to intense bursts of star formation. Yet, the av- 
erage enhancement in star formation rate (SFR) in major mergers and the contribu- 
tion of such events to the cosmic SFR are not well estimated. Here we use pho- 
tometric redshifts, stellar masses and UV SFRs from COMBO-17, 24//m SFRs from 
Spitzer and morphologies from two deep Hubble Space Telescope (HST) cosmological 
survey fields (ECDFS/GEMS and A901/STAGES) to study the enhancement in SFR 
as a function of projected galaxy separation. We apply two-point projected correla- 
tion function techniques, which we augment with morphologically- selected very close 
pairs (separation < 2") and merger remnants from the HST imaging. Our analysis con- 
firms that the most intensely star-forming systems are indeed interacting or merging. 
Yet, for massive (M* > IO'^Mq) star- forming galaxies at 0.4 < z < 0.8, we find that 
the SFRs of galaxies undergoing a major interaction (mass ratios < 1 : 4 and separa- 
tions < 40 kpc) are only 1.80 ±0.30 times higher than the SFRs of non-interacting 
galaxies when averaged over all interactions and all stages of the interaction, in good 
agreement with other observational works. Our results also agree with hydrodynami- 
cal simulations of galaxy interactions, which produce some mergers with large bursts 
of star formation on ~ 100 Myr timescales, but only a modest SFR enhancement when 
averaged over the entire merger timescale. We demonstrate that these results imply 
that only < 10% of star formation at 0.4 < z < 0.8 is triggered directly by major merg- 
ers and interactions; these events are not important factors in the build-up of stellar 
mass since z= 1. 

Subject headings: galaxies: general — galaxies: evolution — galaxies: star formation 
— galaxies: interactions — infrared: galaxies 



Introduction 



Observational evidence from a variety of angles indicates t hat galaxy interactions and merg 



ers of galaxies can lead to dramatically-enhan ced star formation dSanders et al.ll 19881 : iBarton et al 



2000l : iLambas et al.ll2003l : iBarton et al.ll2007|) . This appears to hold true at all redshifts where one 
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can r ecognize mergers 



gies; 



Melbourne et al. 2005: Hammer et al 



through galaxy morpholo gies (z < 1 with rest-frame optical morpholo- 

< z < 3 using less certain UV morpholo- 



rpnoiogie 
I boOSI : 1 



gies lChapman et alj|2004l) . Ultra-luminous infrared galaxies (ULIRGs), representing the highest 



intensity star format ion events at low redshifts, are almost invariably hosted by merging galaxies 
dSanders et al.ll 19881) . For a number of applications the quantity of interest is the average enhance- 
ment in star formation (SF) triggered by merging (ensemble average over the population of major 
mergers/interactions, or equivalently, tem poral average over major merger events during a merger 



temp 



lifetime), not the high-intensity t ail (e.g.. IBarton et al.l 120001 : iLambas et al.l 120031 : iLin et al.l 12007 



Li et al.ll2008] : Ijogee et allboool) . barton et all hooi) carefully quantified the star formation rate 
(SFR) enhancement in mergers in low-mass halos at low redshift, using the Two-Degree Field 
Galaxy Redshift Survey (|Colless et al.ll200lh . They found that roughly 1/4 of galaxies in close 
pairs (separated by < 50kpc) in low-mass halos with M^, < -19 have SFR enhancements of a 
factor of five or moreQ. 

It has also been noted that the strong decrease of the cosmic SFR density betwe en z = 1 
and z = (e.g., IlIUv et al.l[l996l : iMadau et al.ll 19961 : lHopkinsll2004l : Ilc Floch et aPboOSI) was not 
dissimilar from the relatively rapid dro p in merger rate infe rred (at that time) from close pairs 
and morphologically- selected mergers (LeFevreetallboOoi) . If much of the star formation at 
z > 0.5 were triggered by merging, the apparent similarity in evolution between SFR and merger 
rate would be a natural consequence. More recently, studies of the fraction of star formation 
in morphologically-selected interacting and merging galaxies at intermediate redshifts z < 1 have 
demonstr ated that, in fact, the bulk of star formation is in quiesciently star- forming disk-dominated 



galaxies (iHammer et al.ll2005l : IWolf et al.ll2005l : iBell et al.ll2005l : iJogee et al.ll2009D . 



3 



Similarly, it has long been argue d that early-type (elliptica l and lenticular) galaxies ar e a natu- 
ral outcome of galaxy mergers (e.g., iToomre & Toomrdl 19721 : ISchweizer & Seitzenl 19921) . In any 
hierarchical cosmogony mergers are expected to play a large role; a wide range of work — obser- 
vations of the increasing number density of non- star-forming e arly-type galaxies from z = 1 to the 
present (|Bell et al.ll2004l : iBrown et al.ll2007l: iFaber et al.ll2007[). the kinema tic and stellar popula- 
tions of local early-type galaxies ( Trager et al. 2000l : Emsellem et aL 2004 ), or the j oint evolution 



of the stellar mass function and s tar formation rates of galaxies (IBell et al.l 120071 : IWalcher et al 
20081 : IPerez-Gonzalez et al.l |2008|) — has given support to the notion that at least some of the 



early-type galaxies assembled at z < 1 have done so through galaxy merging. In such a picture, 
the average SFR enhancement from merging is of interest for interpreting the SF and chemical 
enrichment history of early -type galaxies, inasmuch as it gives an idea of what kind of fraction of 
stars in present-day early-type galaxies we can expect to have formed in the burst mode, and what 



'This correspo nds rough l y to a m ass cut of 5 x lO^M©, assuming a stellar M/L^j ^ 1, appropriate for a star-forming 
blue galaxy with a lchabriej hoosi) stellar IMF. 
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fraction we can expect to have formed in a quiescent mode in the progenitor galaxies. 



Direct observational constraints on the enhancement in SFR caused by merging provide an 
important calibration for modeling triggered star formation in cosmologically-motivated galaxy 
formation models. Hydrodynamic simulations of interacting galaxies in which gas and star forma- 
tion are explicitly modeled have demonstrated that torques resulting from the merger can efficiently 
strip gas of its angular mor nentum, driving it to high densities and leading to significant enhance- 
ment in star formation (e.g. jBames & Hemquistl[l996l : iMihos & Hemquisil996l : ICox et al.ll2006a[ 
2008l : ldi Matteo et al.ll2007|) . However, state-of-the-art cosmological simulations lack the dynamic 
range to accurately simulate the internal structure of galaxies in significant volumes, so estimates 
of the global implications of merger-driven star formation enhancement have had to rely on semi- 
analytic calculations (e.g. ISomerville et al.ll200ll : iBaugh et al.ll2005l : ISomerville et al.ll2008h . Fur- 
thermore, as the progenitor properties play a key role in the simulated SFR enhancements (e.g.. 



di Matteo et al. 2007: Cox et al. 



20081), inaccurate progenitor property values (e.g., incorrect gas 



fraction or internal structure) will lead to incorrect estimates for the average fraction of SF in 
mergers even if the SF in each individual merger were modeled perfectly. 

Therefore, to constrain galaxy evolution models and to understand the physical processes re- 
sponsible for the main mode of star formation at z < 1 , it is of interest to determine observationally 
the typical enhancemenj^ in SFR averaged over the duration of the entire major (stellar-mass ratio 
between 1:1 and 1:4) galaxy merger or interaction and to constrain the overall fraction of SIF trig- 
gered by mergers/interactions at intermediate redshift. In a companion paper Jo ge e et aL 2009) 
we focus on the rate of merging and also present a preliminary exploration of the average change 
in the SFR caused by late-stage major and minor merging, finding an average mild enhancement 
within the restrictions imposed by the sample size. In this paper we present a statistically-robust 
analysis of the properties of star-forming galaxies at 0.4 < z < 0.8 including all relevant merger 
phases and aimed at providing a satisfactory answer to two key questions. What is the average 
enhancement in star formation rate as a function of galaxy pair separation compared to their SFR 
before the interaction? What fraction of star formation is directly triggered by major mergers and 
interactions? 

There are a number of conceptual and practical challenges in such an experiment. Enhance- 
ments in SFR produce both a boost in luminosity, but also increase dust content and extinction. 
At a minimum, one therefore needs dust-insensitive SFR indicators. In addition, simulations have 
indicated that SF can be e nhanced at almost all phases of an interaction from first passage through 
to after coalescence (e.g.. iBames & Hemquistill996l : Idi Matteo et al.ll2007|) : although close pairs 



^When we refer to SFR enhancement, we define this as the ratio of SFR in some subsample (e.g., close pairs) to 
the average SFR of all systems in that mass bin. 
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will inevitably include some fraction of galaxies before first pass and galaxies with unbound or- 
bits. Therefore, an analysis needs to include both close pairs of galaxies (those before coalescence) 
and morphologically-classified mergers (primarily those near or after coalescence). Morpholog- 
ical classification is not a straightforward art (see ^Jog ee et al. 2009 . for a compariso n between 
automated classifications and visual morphologies), even in ideal cases (lLiskeiil2008|) . Finally, 
galaxy mergers are rare and short-lived, necessitating large surveys to yield substantive samples of 
mergers. 



In this w ork, we address these challenges as far as possible (see also Lin et al.l 120071 and 



Li et al.l2008h. We use e stimates of redshift and stellar mass from the COMBO- 17 survey (IWolf et al 



20031 : iBorch et al.ll2006|) to define and characterize the sample. Stellar mass selection should limit 



the effect of enhanced star formation and dust content on the sample definition. We use SFR indi- 
cators that are constructed to be dust extinction insensitive, by combining ultraviolet (UV; direct, 
unobscured light from young stars) and infrared (IR; thermal emissio n from heated d ust, pow- 
ered primarily by absorption of UV light from young stars) radiation (IBell et al.ll2005|) . Finally, 
we study a very well-characterized sar nple of galaxy pairs at 0.4 < z < 0.8 using weighted pro- 
jected two-point correlation functions (|Skibba et al.ll2006l : iLi et al.ll2008h . supplementing them at 
very small separations < 15kpc with very c lose pairs or rn erger remnants morphologically se- 
lected from two wide HST mosaics, GEMS (IrIx et all 120041) and STAGES (Gray et allboooh . in 
an attempt to account for all stages of galaxy interactions. 

The plan of this paper is as follows. In ^we discuss the data and the methods used to estimate 
the stellar masses and the SFRs. In §|3]we describe the sample selection and the method used for 
the analysis. In ^we present our estimates of the enhancement in SFR as a function of projected 
separation. In §|5J we compare with previous observations, constrain the fraction of SF triggered by 
major mergers and interactions at 0.4 < z < 0.8, and compare with simulations of galaxy merging. 
Finally in ^we summarize the main findings of this paper. All the projected distances between the 
pairs used here are proper distances. We assume Hq = 70kms"^ Mpc"\ Qao = 0-7 and = 0-3. 



2. The Data 



2.1. COMBO-17. Redshifts and stellar masses 



COMBO- 17 has to date fully surveyed and analyzed three fields to deep limits in 5 broad 
and 12 medium p ass- bands (Extended Ch andra Deep Field South (ECDFS), A901/2 and SI 1, see 
Wolf et al.l (|2003|) and lBorch et al.l (|2006|) ). Using galaxy, star and quasar template spec tra, objects 
are classified and redshifts assigned for ~ 99% of the objects to a limit of nin ~ 23.5 (|Wolf et al 
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2004|) . The photometric redshift errors can be described as 



l+Z 



0.007 X [l + lo°•*^(""^-^^■''Y/^ 



(1) 



and rest frame colors and absolute magnitudes are accurate to ~ 0. 1 mag (accounting for distance 
and ^-correction uncertainties). The astrometry is accurate to ~ 0. l"and the average seeing is 0.7". 
It is worth noting that Eq.[T]leads to typical redshift errors of ~ 0.01 for bright (m« < 21) and 
~ 0.04 for faint (21 < < 23.5) galaxies in the 0.4 < z < 0.8 interval. 



The stellar masses were estimated in COMBO- 17 by lBorch et al.l(|2006h using the 17-passband 
photometry in conjunctio n with a non-evolving template library derive d using the Pegase stel- 
lar population model (see lFioc & Rocca-Volmerangd 1 19971119991) and alKroupa et al.l (Il993b ini- 



Chabriei 



( 20031) 



tial mass function (IMF). Note that the results assuming a iKroupal (|2001|) or a 
IMF yield similar stellar masses to within ~ 10%. The reddest templates have smoothly-varying 
exponentially-declining star formation episodes, intermediate templates have a contribution from 
a low-level constant level of star formation, while the bluer templates have a recent burst of star 
formation superimposed. 

The masse s are consistent with those using M/L estimates based on a single color (e.g.. 
Bell et al.l 120031) . Random stellar mass errors are < 0.3 dex on a galaxy-by-galax y basis, and 
systematic errors i n the stellar masses we re argued to be at the 0. 1 dex level (see B orch et al.ll2006L 
for more details). iBelL&^de JongI feoOll) argued that galaxies with large bursts of recent star for- 
mation could produce stellar M/L values at a given color that are lower by up to 0.5 dex; this 
uncertainty is more relevant in this work than is often the case. While this will inevitably remain 
an uncertainty here, we note that the lBorch et al.l(|2006|) templates do include bursts explicitly, thus 
compensating for the worst of the uncertainties introduced by bursting star formation histories. In 
§ 14. 1.21 we will explicitly study the impact that such uncertainties have on our results. 

In what follows, we use COMBO- 17 data for two fields: the ECDFS and Abell 901/902 fields, 
because of their complementary data: deep HST/ACS imaging from the GEMS and STAGES 
projects respectively (allowing an investigation of morphologically-selected merger remnants and 
very close pairs), and deep 24/im imaging from the MIPS instruments on board Spitzer, required 
to measure obscured SF. 



2.2. GEMS and STAGES HST imaging data 



F606W (V-band) imaging from the GEMS and STAGES surveys provides 0.1" resolution im- 
ages for our sam ple of COMBO- 17 galaxies. Using the Advanced Camera for Surveys (ACS; 
Ford et al.ll2003b on board the Hubble Space Telescope (HST), areas of ~ 30' x 30' in each of the 
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ECDFS and the A90 1/902 field have bee n surveyed to a depth allowing galaxy detection to a lim - 
iting magnitude of mfi^,{F606W) = 28.5 JrIx et alJl2004l : lOrav et alJl2009l : Icaldwell et aUboOsI) . 
These imaging data are later used to visually classify galaxies, allowing very close pairs (separa- 
tions < 2") and merger remnants to be included in this analysis. We choose not to use F850LP 
HST data available for the GEMS survey in order to be consistent in our classification between the 
two fields (only F606W is available from STAGES). 



2.3. MIPS 24 /ini , total infrared emission and star formation rates 



The IR observatory Spitzer has surveyed two of the COMBO-17 fields: a 1° x 0.°5 scan of the 
ECDFS (MIPS GTO), and a similarly-sized field around the Abell 901/902 galaxy cluster (MIPS 
GO-3294: PI Bell). The final images have a pixel scale of 1^^25/pi xel and an image PSF FWHM of 
~ 6". Source dete ction and p hotometry are described in depth in IPapovich et al.l (120041) and cata- 
logue matching in Bell et al. ( 2007^ Based on those works, we estimate that our source detection 
is 80% complete at the 5a limit of 83/iJy in the 24yum data in the ECDFS for a total exposure 
of ~ 1400 s pix"^ The A90 1/902 field has similar exposure time, but owing to higher (primarily 
zodiacal) background the 5a limit (80% completeness) is 97yuJy, with lower completeness of 50% 
at 83^Jy. We use both catalogs to a limit of 83/iJy. 

To include both obscured and unobscured star formation into the estimate of the SFR of 
galaxies in our sample, we combine UV emission with an estimate of the total IR luminosity 
in concert. As the total thermal IR flux in the 8-1000/im range is observationally inaccessible 
for almost all galaxies in our sample, we have instead estimated total IR luminosity from the 
observed 24/im flux, corresponding to rest-frame 13-17/im er nission at the redshifts of interest z = 
0.4-0.8. For this exercise, we adop t a Sbc template from the iDevriendt et al.l (|1999h SED library 
dzheng et al.l 120071 : iBell et al.ll2007b . The resulting IR luminosity is accurate to a factor of < 2: 
local galaxies with IR luminosities in excess of > l O^^Lp, show a tight co r relation between rest- 
frame 12-15/im luminosity and total IR lu minosity dSpinoglio et al.lll995l : ICharv & Elbazll2001.: 



Roussel et al. l l200ll : IPapovich & Beli2002h with a scatter of ~ 0. 15 dex. Furthermore, Zheng et al . 
J2OO7I) have stacked luminous (Ltir > IO^^Lq) z ~ 0.7 galaxies at 70^m and 160Atm , find i ng tha t 



their average spectrum is in good agreement with the Sbc template from IDevriendt et al.l (|l999h . 
validating at least on average our choice of IR SED used for extrapolation of the total IR luminosity. 



We estimate the SFR by using both directly observed UV-light from massive stars and dust- 



bin this paper, we are interested in SFR enhancements in close pairs of galaxies, where the closest pairs may fall 
within a single 5f>iteer/MIPS PSF. Accordingly, in this work we choose to explore the total SFR in the pair (avoids 
deblending uncertainties) rather than the individual SFR occuring in both galaxies. 
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obscured UV-light measured from the mid-infrared. As in lBell et alj (|2005|) we estimate the SFR tp 
by means of a calibration derived fr om Pegase syn thetic models assuming a 100 Myr-old stellar 
population with constant SFR and a lChabrieii (|2003|) IMF: 



V'/CMeyr-i) = 9.8 x 10"" x (Ltir + 2.2Luv). (2) 

Here Ltir is the total IR luminosity and Luv = l-5z//,^ 28oo is a rough estimate of the total in- 
tegrated 1216A-3000A UV luminosity. This UV luminosity has been derived from the 2800A 
rest-frame luminosity from COMBO- 17 /y,28oo- The factor of 1.5 in the 2800A-to-total UV con- 
version accounts for UV spectral shape of a 100 Myr-old population with constant SFR, and the 
UV flux is multiplied by 2.2 to account for the light emitted longwards of 3000A and shortwards 
of 1216A by the unobscured stars belonging to the young population. 

For all galaxies detected above the 83/xJy limit, we have used the IR and UV to estimate the 
total SFR. For galaxies undetected at 24//m, or detected at less than 83//Jy, we use instead UV-only 
SFR estimates. 



2.3.1. IR emission from AGN -heated dust 



Possible contamination of mid-IR-derived SFRs from AGN heated dust is often addressed by 
estimating the fraction of star formation held in X-ray detected sources. In our case < 15% of the 
sta r forming galaxy sample were dete cte d in X-rays , in go od agreement with the results found by 



I.e. 



Silva. Maiolino. & Granatol (120041) or lBell et all (120051) 



Yet, there are two limitations of this estimate. Firstly, this does not account for any con- 
tribution from X-ray undetected Compton-thick AGN, which could drive up the expected con- 
tribution from AGN in our sample. For example, applying an mn = 24 cut to the sample of 
(l2006l) . we estimate the fraction of X-ray undetected AGN to be ~ 30%, 
(119991) find ~ 50% of local AGN to be Compton thick. On this basis, it is 



Alonso-Herrero et al 



while 



Risaliti et al. 



conceivable that up to 30% of 24/im luminosity is from galaxies with AGl^ , 

Secondly, even in galaxies with AGN, not all of the IR emission will come from the AGN. 
Although the data does not currently exists to answer this question conclusively, it is possible 
to make a rough estimate of the effect. In order to estimate the fraction of mid-IR light that 



^Alth ough note that in a recent investigation of X-ray undetected IR-bright galaxies in the CDFS, iLehmer et al. 
din prep J) found that radio-derived (1.4GHz) SFRs agree with the UV+IR-derived ones. This implies that the relative 
strength of any AGN component is not dominant when compared to the host galaxy. 
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co mes from the AGN (as oppos ed to star formation in the host), we have made use of the results 



of iRamos Almeida et alj (|2007|) . who attempted to structurally decompose mid-infrared imaging 
from Infrared Space Observatory for a sample of both Seyfert 1 and 2 AGN in the local Universe, 
some of which are very hig hly-obscured in X-rays. Analyzing the results in Tables 2 and 3 of 



Ramos Almeida et al.l (|2007[) . we have found that only a small fraction of the IR radiation at ~ 
10/im (in this paper we work at rest-frame 13-17/im ) comes from the central parts of the galaxies 
in the Seyfert 2 population, finding a total contribution of: 



?AGN 

=0.26 ±0.02. (3) 



fptotal 
^IR 



This result should be viewed as indicative only: obviously, the systems being studied will be 
different in detail from those in our sample. Furthermore, the lO^m luminosities of the nuclei will 
be preferentially affected by silicate absorption, making it possible that our value of Fjr^ /Fj'^'"' is 
a lower limit. 

Despite the various levels of uncertainty, taking the different lines of evidence together demon- 
strates that < 30% of the IR luminosity in our sample comes from systems that may host an AGN, 
and that it is likely that < 10% of the IR luminosity of our sample is powered by accretion onto 
supermassive black holes. Given the other uncertainties in our analysis, we choose to neglect this 
source of error in what follows. 



3. Sample selection and method 

The goal of this paper is to explore the star formation rate in major mergers between massive 
galaxies, from the pre-merger interaction to after the coalescence of the nuclei. We chose a stellar 
mass-limited sample with > lO^^M© in the redshift slice of 0.4 < z < 0.8 (see Fig. 1). This 
roughly corresponds to My = -18.7 for galaxies in the red sequence and My = -20.1 for blue 
objects. We only included galaxies that fall into the footprint of both the ACS surveys GEMS and 
STAGES and of existing Spitzer data. These criteria resulted in a final sample of 2551 galaxies. 

Given the fl ux limit mr < 23.5 for which COMBO- 17 has reasonably complete redshifts 
(|Wolf et al.l 120041) we are complete for > IO^^A/q blue cloud galaxies over the entire redshift 



range 0.4 < z < 0.8. For red sequence galaxies, the sample becomes somewhat incomplete at 
z > 0.6, and at z = 0.8, the limit is closer to 2 x 1O"'M0. We chose to adopt a limit of IO^^Mq 
in what follows, despite some mild incompleteness in the red sequence, for two reasons. First, 
adopting a cut of 2 x IO^'Mq across the whole redshift range reduces the sample size by a factor of 
30%, leaving too small a sample for the proposed experiment. Second, the vast majority of the star 
forming galaxies are blue cloud galaxies (83% of the star formation is occuring in blue galaxies), 
making the modest incompleteness in the red sequence of minor importance. 
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Later, we will use a subsample composed of star forming galaxies. We will refer to 'star 
formers' as galaxies defined by having either blue optical colors or having been detected in the 
MIPS 24//m band. We select opticaUy-blue ga laxies adopting a stellar mass-dependent cut in rest- 
frame U-V color, following Bell et all J2OO7I) : U -V > 1.06-0.352z+0.227(logioM,- lofl (see 
Fig.[T]). We include all objects detected above the lAfim limit of 83/i Jy as star-forming. 

In order to track star formation in very close pairs (< 2" and hence unresolved by the ground- 
based COMBO- 17 data) and merger remnants, we include only merging systems (from the ACS 
data) with > 2 x 10"'Mq: i.e. the minimum possible mass for a merger between two galaxies 
in our sample. 




Fig. 1. — Stellar mass vs. color distribution of COMBO- 17 selected galaxies in the ECDFS and 
A901/2 field with 0.4 < z < 0.8. The vertical line shows the mass limit > IO'^Mq used to select 
our sample. This mass selected sample is complete except for red sequence galaxies at z > 0.6. 
The blue line shows the cut used to separate red sequence and blue cloud galaxies. Red symbols 
denote 24 fim detected galaxies with >S3fi Jy. 



^Due to minor magnitude and color calibration differences between the two fields, the red sequence cut is slightly 
field dependent, with the intercept at IO^'Mq and z=0 being U-V=1.01 and 1.06 for the ECDFS and the A901/902 
fields. 
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3.1. Projected correlation function 



The correlation function formalism is a conven i ent and powerful tool to c haracterize pop- 
ulations of galaxy pairs (e.g. bavis & Peebleslll983l : lieisbart & Kerscheilbooci) . Here, we use 
weighted projected two-point correlation functions because redshift uncertainties (1-3%) from 
COMBO- 17 translate to line-of-sight distance errors of ~ 100 Mpc, necessitating the use of pro 



jected correlation functions to explore the properties of close physical pairs of galaxies (IBell et al. 



20061). For our sarnple at hand, we estimate the weighted (or marked) two-point correlation func 



tion (Boerner et al 



19891 : iBeisbart & Kerscheil l2000l : ISkibba et al.ll2006l : ISkibba & ShethI 120081) . 



using both the SFR and the specific SFR (SFR per unit stellar mass) as the weight. 

The projected correlation function w(rp) is the integral along the line of sight of the real-space 
correlation function: 



wirp)= / a[rl + 'Jr^]'hdn, (4) 



where is the distance between the two galaxies projected on the plane of sky and tt the line-of- 
sight separation. A simple estimator for this unweighted correlation function is w{rp) = A(DD / RR- 
1), where A is the path length being integrated over, DD{rp) is the histogram of separations be- 
tween real galaxies and RR(rp) is the histogram of separ ations between g alaxies in a randomly- 



distributed catalogue (this is the same estimator used in iBell et al.l 120061 ). Basically, the aim is 
to find the excess probability (compared to a random distribution) of finding a galaxy at a given 
distance of another galaxy. This estimator accomplishes that by subtracting the random probabil- 
ity of finding two galaxies at a given separation from the probability in the real data sample and 
normalizing to the probability in the random case. Other estimators (i.e. A[{DD - DR) / RR] or 
A[(DD-2DR+RR)/RR]) for the 2-point correlation function give results different by < 5% (less 
than other sources of uncertainty). Thus: 



DD(rp) = J2Di 

ij 

RR(rp) = J2Rij 



'J 

where the sum is over all non-repeated pairs in the sample, and D,y (Rij) equals 1 only if the pair 
selection criteria are satisfied in the real (random) galaxy catalogue, and is equal to otherwise. 
The first criterion is that the stellar-mass ratio falls between 1:1 and 1:4. We further only allow 
a maximum redshift difference Az = A = V2a,, where cTj is the error in redshift of the primary 
galaxy (see Equation [B, and, depending on the case, either the primary or both galaxies in the pair 
have to be star formers (see 
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We can then study the possible enhancement of (specific) star formation rate by means of a 
projected marked (or weighted) correlation function, which can be defined: 



1 + W^(r„)/A 
l+w(rp)/A 



where W{rp) = A{PP/PPr- 1) and, 



PP(rp) = Y/ipij 

ij 

PPRirp) = Y/uRij- 

u 

Pij is the mark (or weight). We adopt two different weights Pij in what follows, one is the SFR of 
the pair of galaxies: 

Pij = = SFRij = SFRi + SFRj, 
and the other is the specific SFR of the galaxy pair: 

o r.r.r, SFRi + SFRj 

P, = ., = Specific5F7?,= ^^^^^^^ ■ 
Then, the estimator that we use for E(rp) is: 

PP/DD 

£W = ^, (6) 

where {Ptj) is the average value of the weight used (SFR, or specific SFR) across the sample. It 
is worth noting that with our definition of the enhancement given in eq. 5 the random histograms 
RR and PPr cancel in the process of obtaining the expression in eq. 6, so they are not used in the 
computation of our enhancement. 

In the present work we perform two analyses: the cross-correlation of star forming galaxies (as 
defined above) as primary galaxies with all galaxies as secondaries, and the autocorrelation of star- 
forming galaxies. We will estimate the errors in our mark by means of bootstrapping resampling. 



3.2. Visual Morphologies 

A particular challenge encountered when constructing a census of star formation in pairs and 
mergers is accounting for systems with separations of < 2" (which corresponds to < 1 5 kpc, the 



radius within which we can no longer separate two massive galaxies using COMBO- 17: iBell et al 



20061). In order to pick up the SF in all the stages of the interaction, we need to have an estimate 
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of the SFR not only in galaxy pairs with separations > 15 kpc but also in extremely close pairs and 
in recent merger remnants. We conduct our census of such close physical pairs by including in the 
< 15 kpc range sources that are not resolved by COMBO-17, but appear to be interacting pairs or 
merger remnants on the basis of visual classification of the ~ 0. 1" resolution ACS images. We try 
to recover visually all < 15 kpc separation pairs of two > IO^^^Mq galaxies with a mass ratio 
between 1:1 and 1:4 missed by COMBO-17. In addition to those extremely close pairs, we also 
account for the SF in recent merger remnants > 2 x 10^°A/q (two times the minimum mass of 
a galaxy in the sample and the minimum possible mass of a galaxy pair as defined before). 



3.2.1. Discussion of visual classifications 



Our goal is to include very close pairs or already-coalesced major merger remnants into the 
census of 'mergers' in order to account for any SF triggered by the merger/interaction procesflwe 
do so on the basis of visual classification of the sample. The motivation for visual classification 
is a pragmatic one: while a number of automated m orpholo gical classification systems have been 
developed in the last 15 years (i.e. Abraham et al. 19961 : Iconselice et al. 2003 : Lotz et al. 2004 , 
etc.), it seems that the sensitivity of the observables used (asymmetry, dumpiness, Gini coefficient, 
second order moment of the 20% brightest pixels) is insufficient for matching the performance of 
visual classification in current intermediate redshift galaxies with th e same level of precision that 



they display in the loca l Universe samples used for their calibration (|Conselice et al 
20081 : Ijogee eTaDl2009l) . 



20031 : iLisker 



Yet, there is a degree of subjectivity to what one deems to be a major merger remnant. Many 
factors shape the morphology of a galaxy merger that are beyond the control of the classifier. 
Bulge-to-total (B/T) mass ratios have an strong effect on both the intensity of the SFR enhance- 
ment and the time at which the intensity peak shall occur (e.g. jMihos & Hernquistill996h . Orbital 
parameters strongly shape the development of easily recognizable tidal tails and bridges (copla- 
nar or not, retrograde vs. prograde, etc). Prior dust and gas content of the parent galaxies ('dry' 
vs. 'wet' mergers) will make a difference to the appearance of the final object during the coa- 
lescence. Furthermore, merging timescales will depend on whether the galaxies are undergoing 
a first passage or are in the final stages of the merger. Finally, there is a degeneracy between all 
these parameters and the relative masses of the galaxies undergoing the interaction, which makes 
difficult in some cases to distinguish the morphological signatures of a major merger from those 



^ Note that a consistent comparison with the projected correlation function sample requires the inclusion of all 
non-interacting pairs that are physically-associated (in the same cluster, filament, etc), are seen to be close projected 
pairs on the sky, but may be separated by as much as a few Mpc along the line of sight. 
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of a minor merger. 



Some of these factors (e.g. gas fraction, B/T ratios, etc) w i ll also affect the enhancement of 
the SFR during the interaction (e.g. Idi Matteo et allboOTL boosi : Icox et allboos l). While there is 
considerable merger-to-merger scatter, encounters of two gas rich disk galaxies with parallel spins 
tend to develop, on average, the strongest morphological features, but at the same time are more 
likely to throw out large amounts of cold gas in tidal tails, preventing the funneling of this gas to 
the central regions. Thus, samples selected to have the strongest morphological features may have 
an average SFR enhancement different from the actual mean enhancement- 
One practical issue is that of passband choice and shifting. We choose to classify the F606W 
images of the GEMS and STAGES fields (in STAGES because that is the only available HST pass- 
band and in GEMS for consistency and because F606W has higher S/N that the F 850LPdata). This 
corresponds to rest-frame ^ 430(3 30)nm at redshift 0.4(0.8). In previous papers ( IWolf et al.ll2005 



Bell etaLl llOOslJjogee et al.l 120091) . we have assessed whether the morphological census derived 
from GEMS/STAGES would change significantly if carried out data a factor of 5 deeper from the 
GOODS project (testing sensitivity to surface brightness limits), or if carried out at F850LP (always 
rest-frame optical at these redshifts). We found that the population does not show significantly 
different morphologies between our (cor nparatively) shallo w F606W data and the deeper/redder 
imaging data from GOODS (see Fig. 5 in ljogee et al.ll2009f) . 



3.2.2. Method 

An independent visual inspection of the galaxy sample has been carried out by four classi- 
fiers, A.R.R., E.F.B., R.E.S. and D.H.M. in order to identify morphological signatures of major 
gravitational interactions. Each classifier assigned every one of the ~2500 sample members to one 
of the three following groups: 

1. Non-major interactions: The bulk of galaxies in this bin show no signatures of gravitational 
interactions. Asymmetric, irregular galaxies with patchy star formation triggered by internal 
processes lie in this category. A small fraction of galaxies in this bin show a clearly recog- 
nizable morphology (e.g. spiral structure) but also signatures of an interaction (such as tidal 
tails, or warped, thick or lopsided disks) but have no clear interaction companion; note that 
these objects could be interacting systems where the companion is now reasonably distant 
and/or faint and more difficult to identify. The tidal enhancement of SF from such systems 



^This bias might also be present in the case of studies looking for signs of interactions in the host galaxies of AGNs, 
attempting to assess whether the AGN activity is preceded by a merger. 
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will not be missed by putting them in this bin; rather, it will be measured statistically and 
robustly from the two point correlation function analysis. Minor mergers and interactions 
(interactions where the secondary is believed, on the basis of luminosity ratio, to be less than 
1/4 of the mass of the primary) also belong to this category. 

2. Major close interactions: Close pairs resolved in HST imaging but not in ground-based 
COMBO- 17 data, consisting of two galaxies with mass ratios between 1:1 and 1:4 based 
on relative luminosity, and clear signatures of tidal interaction such as tidal tails, bridges or 
common envelope (see Fig.[2l). From now on we shall refer to objects classified in this group 
as "very close pairs". 

3. Major merger remnants: Objects that are believed to be the coalesced product of a recent 
major merger between two individual galaxies. Signposts of major merger remnants include 
a highly-disturbed 'train wreck' morphology, double nuclei of similar luminosity, tidal tails 
of similar length, or spheroidal remnants with large-scale tidal debris (see Fig. [3]). Galaxies 
with clear signs of past merging but a prominent disk (e.g., highly asymmetric spiral arms or 
one tidal tail) were deemed to be minor merger remnants and were assigned into the group 
1. Naturally there is some uncertainty and subjectivity in the assignment of this class, in 
particular; such uncertainty is taken into account in our analysis by the monte Carlo sampling 
of all four classifications in order to properly estimate the dispersion in the opinions of the 
individual classifiers (see below). 

Table 1 : Results from the morphological classification 



Lower mass limit 


Sample size 


Group 1 


Group 2 


Group 3 


10"'Mq 


2551 


2380 ±37 ±49 


106±7±10 


72±7±8 


2 X 1O^°M0 


1749 


1640 ±32 ±40 


69±6±8 


44±5±7 



Note. — Galaxy and interaction sample. Group 1: Isolated objects and minor interactions. Group 2: Extremely 
close pairs (rp < 15 kpc). Group 3: Merger remnants. The first error bar represents classifier-to-classifier scatter while 
the second one represents Poisson noise. 

We then assign the objects in the groups 2 and 3 (very close pairs with morphological sig- 
natures of interaction and merger remnants respectively) to a small projected separation and treat 
every one of them as a galaxy pair in order to combine them with the correlation function analysis 
result for pairs with separations > 2". All objects in Group 2 (extremely close pairs with projected 
separations < 15 kpc as measured by centroids in HST imaging) are assigned to a separation of 
10 kpc and all objects in Group 3 (merger remnants) are assigned to a separation of kpc. We 
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have checked for duplicate pairs in both the visually selected sample and the COMBO- 17 catalog 
in order to avoid repeated pairs. Galaxies in group 1 are already included in the two point corre- 
lation function analysis, and any SF triggered by major interactions or early-stage major merging 
is accounted for by that method. As we have four different classifications for every object (one 
given by each human classifier), we randomly assign one of them, calculate the average value of 
the weight we are using and repeat the process a number of times. This approach presents two 
clear advantages: a) the resultant bootstrapping error not only represents the statistical dispersion 
but also the different criteria of the four human classifiers, and b) the morphology of every object is 
weighted with the four classifications given. This means that objects with discrepant classifications 
are not just assigned to one category when we calculate the SFR (or specific SFR) enhancement; 
rather, any dispersion in classifications is naturally accounted for (e.g. minor/major criteria). The 
numbers of such systems and their uncertainties, estimated from the classifier-to-classifier scatter, 
are given in Table [T] 



4. Results 

We are now in a position to quantify the triggering of star formation in galaxy interactions and 
mergers in the redshift interval 0.4 < z < 0.8, in the cases where each galaxy has > IO^^Mq and 
the pair has a stellar mass ratio between 1:1 and 1:4. Our primary analysis is based on a marked 
cross -correlation between star-forming galaxies, as defined in §3, and all galaxies in the sample. 
For morphologically-selected very close pairs or interactions (unresolved by COMBO- 17), we also 
require them to be blue or detected by Spitzer to be considered as part of the star- forming sample^, 
with a mass of > 2 x IO'^Mq. 

We perform two analyses in this paper: the cross -correlation of star-formers as primary galax- 
ies with all galaxies as secondaries (our default case), and the autocorrelation of star-forming 
galaxies. While the first analysis is a rather more direct attack on the question of interest, we show 
results from the autocorrelation of star- forming galaxies to illustrate the effects of making different 
sample choices on the final results. 




Fig. 2. — Objects classified in group 2: Major close interactions. The presence of two galaxies 
and signs of interaction are required. The classifier believes the mass ratio is between 1:1 and 
1:4. At this stage of the interaction, dry mergers are still recognizable as seen in panels at top 
center, bottom center and bottom right. The black bar at the bottom of every panel shows a proper 
distance of 20 kpc at the redshift of the object. Some of the objects classified in this group were 
also separated as two galaxies in the ground-based catalog and treated in consequence. 
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Fig. 4. — Pair specific SFR enhancement as function of the projected separation between two 
galaxies. The two smallest radii bins are derived from morphologically-selected very close pairs 
(shown with rp ~ 10 kpc) and merger remnants (shown with rp = 0); enhancements at larger radii 
are determined using weighted two-point correlation functions. A statistically significant enhance- 
ment is present in galaxy pairs and mergers below 40 kpc in both the cross-correlation between 
star forming galaxies as primaries and all galaxies as secondaries (black filled symbols) and the 
autocorrelation of star-forming galaxies (empty diamonds). Error bars have been calculated by 
bootstrapping. 



4.1. Enhancement in the Star Formation Activity 

Our main results are shown in Fig. |4l which shows the enhancement of the specific star for- 
mation rate (SSFR) in pairs as a function of their projected separation. As explained in ^2.31 we 
use UV+IR SFRs for the objects detected in 24 /^m and only UV SFRs for those undetected. For 
the whole sample, 38% of the galaxies where detected by Spitzer above the 83/xm limit, while if 
we restrict to the groups 2 and 3 in our morphological classification we find a detected fraction of 
60%. Fig. in shows a clear enhancement in the SSFR for projected pair separations rp < 40 kpc. It 



^All galaxies, irrespective of their color or IR flux, were classified; the star-forming galaxies are simply a subsample 
of this larger sample. 



-20- 



could be argued that the SSFR is a better measure of the SF enhancement than the SFR-weighted 
estimator, because the strong scaling of SFR with galaxy mass is factored out. The figure shows 
both the cross-correlation between star-forming primaries and all secondaries (SF-AU, solid line) 
and the star-forming galaxy autocorrelation (SF-SF, dotted line). The two bins at Vp < 15kpc 
are calculated from morphologically-selected very close pairs (r^ = lOkpc) and merger remnants 
(rp = 0). All the errors in Essfr have been computed by bootstrap resampling. This approach al- 
lows us to treat both the morphologically-selected objects and the galaxy pairs exactly in the same 
way, having as a result a coherent display of the error bars. 

There are two reasons why this excess in Essfr in close pairs and remnants is likely a sign 
that interactions induce additional star formation, rather than being due to a correlat ion with 



some other unidentified quantity: a) It is well known from simulations (|Mihos & Hemqu ist 1996; 



di Matteo et al.l 120071: ICox et al.ll2008|) that a burst of star formation is expected in the collisions 



of gas-rich g alaxies, and b) the o bserved effect is in the opposite sense of the usual SFR-density 



relation (e.g.. lBalogh et al.ll2002|) . which says that gal axies in dense envir onments (where preferen 



tially close galaxy pairs tend to be found, as shown in barton et al. I I2OO7I) have, on average, weaker 



star formation activity than galaxies in less dense regions. 

Even when we consider our morphological classification and further Monte Carlo resampling 
method to be very robust, potential classification errors could act in two different directions. Inter- 
acting systems misidentified as non-interacting will be diluted into the background star formation 
as single galaxies contributing to pairs at random separations. While this SF should be lost to the 
interacting bin, the effect on the average SFR would be minimal. On the other hand, isolated galax- 
ies misidentified as interacting systems because of internal instabilities or stochastic star formation 
would act to reduce the enhancement. 

As mentioned before, the SFR for the objects undetected at 24 //m has been calculated based 
only on the UV. In the 24/im detected objects, we have found no clear trend in both the UV 
vs. UV-i-TIR SFRs and in the TIR/TUV vs. optical dust attenuation but found instead a constant 
correction factor with a large scatter (4. 1 ± 2.4 as estimated from the relation between TIR/TUV vs 
optical attenuation.) We have checked the effects of such a dust-correction of the UV-only SFRs: 
the results differ in all bins by <10%, comparable to or smaller than other sources of systematic 
uncertainty. 

Yet, in order to understand the degree of obscuration in galaxy interactions we have repeated 
our analysis including only UV-derived SFRs, this is, excluding the TIR component in Eq. |2]for 
24/im detections. The result of this analysis is shown in Fig. [51 The enhancement in the unob- 
scured SSFR measured for close pairs (rp < 40kpc) in this case is dramatically smaller than the 
enhancement including the dust-obscured (IR-derived) star formation rate. This is more appar- 
ent in the very close pairs and merger remnants, where the excess in the SSFR even disappears 
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completely in the case of the SF-SF autocorrelation {E{rp < 15kpc) ~ 1). This implies that most 
of th e directly triggered s t ar for r nation is dust obscured , in go o d agreement with the expectations 
from lMihos & Hemauistl (Il994[ Il996l) . Idi Matteo etaP toO'^ . Icox et all dlOOsb and the detailed 
models by IJonsson et al.1 (|2006|) . In these simulations most of the star formation is triggered in the 
central regions of the galaxy after the cold gas has been funneled to the inner kpc. 
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Fig. 5. — Same as Fig. ^but tracing only unobscured (UV-derived) star formation. The un- 
obscured SSFR enhancement found in galaxy pairs with separations rp < 40 kpc and mergers 
remnants is dramatically reduced with respect to the case in which the obscured star formation is 
taken into account (Fig. I?])- 



This scenario is also supported by our measurement of the mean ratio between the total SFR 
and the UV-derived, which gives an idea of the degree of dust-obscuration (SFRjp+uv / SFRuv)- We 
find 6.64 ± 0.66 in the case of the merger remnants (Group 3 in ^3.21) and 6.63 ± 0.64 in the case 
of the very close pairs (group 2), compared to 3. 15 ± 0.53 for all objects in the sample. 



4.1.1. Star formation rate vs. specific star formation rate 

To study the fraction of the global star formation directly triggered by galaxy-galaxy interac- 
tions the enhancement in the SFR (rather than in the SSFR) is a better quantity to consider. 

We show in Fig. [6]the enhancement in the SFR (EsFRirp)) as a function of the projected pair 
separation. For the cross-correlation function (our default case) the enhancement in the SFR is 
similar to the one found in the SSFR at all separations except for the merger remnants (r^ = 0), in 



-22- 



which the excess above the whole population is ~ 50% lower. The SFR-weighted autocorrelation 
of star forming galaxies matches that of the SSFR-weighted one for Vp < 40 kpc but differs beyond: 
EsFR = 1 -25 for 40 < rp < 1 80 kpc. While most of these points in Fig.[6]are individually compatible 
with the error bars shown in Fig. IH taken together they represent a ~ 2a significant difference 
between Essfr and Esfr for the entire region 40 < rp < 180 kpc. 




Fig. 6. — SFR enhancement in galaxy interactions. The two smallest radii bins are derived 
from morphologically-selected extremely close pairs (rp ^ \Okpc) and merger remnants (shown 
with rp = 0); enhancements at larger radii are determined using weighted two-point correlation 
functions. There is a clear enhancement at rp < 40 kpc for the cross-correlation analysis (black- 
filled symbols) which is compatible with Essfr (Fig- 0) except for the merger remnants, where the 
excess is ^50% lower The autocorrelation of star forming galaxies (empty symbols) presents an 
unexpected behavior, showing a very mild enhancement at rp < ISO kpc. 



A potential driver of the SFR enhancement in the regime 40 < rp < 18 kpc is the fact tha t 
more massive galaxies tend to be both more clustered and have higher SFR (INoeske et al.ll2007b : 
this could translate into a wea k enhancement in the SFR in galaxy pairs living in dense environ- 
ments (see lBarton et al.ll2007l for a thorough discussion on the relation between galaxy pairs and 
environment) which will not be present in the SSFR, because the normalization by galaxy mass 
factors out this dependence. To test the relevance of this systematic effect, we randomized the 
SFRs among galaxies of similar mass 500 times in the sample and repeated the analysis. We show 
the results of this exercise in Fig. Ul where we can see a tail of enhancement with a behavior similar 
to the one seen in Fig. [6l We believe that a combination of the density-mass-SFR relation plus 
noise is driving Esfr > 1 (autocorrelation) between 40 and 180 kpc. 
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Fig. 7. — SFR enhancement measured after randomizing the SFR between galaxies of similar 
stellar mass. A mild enhancement is found out to separations of ~ 160 kpc. We show the points 
corresponding to the SF-SF autocorrelation at distances > 15 kpc, where no morphological infor- 
mation is used. 

Accordingly, we consider only the enhancement at rp < 40 kpc as produced by major merging 
in what follows, and use the differences between the SFR and SSFR enhancement on < 40 kpc 
scales as a measure of systematic uncertainty. Under those assumptions, we find a weak en- 
hancement of star formation at rp < 40 kpc of e = 1.50 ±0.25 in the SF-SF autocorrelation and 
e = 1.80 ±0.30 in the SF-AU cross -correlation. These values have been computed as the average 
of the enhancement in the bins rp < 40 kpc together in Essfr and Esfr- These (conservative) error 
bars include both the statistical uncertainties and the systematics driving the differences between 



As we have briefly mentioned in ^ there are some uncertainties which need to be estimated 
in the process of calculating the enhancement in the SF activity. Here we try to estimate the impact 
of the stellar-mass and IR SED selection uncertainties. Through this section we will focus in our 
default case, the SF-AU cross-correlation. 



the SFR and the SSFR. 



4.1.2. Further Uncertainties 



Random errors i n stellar masses in iBorch et al.l (|2006|) are < 0.1 dex (with 0.3 dex in cases 
with large starbursts (|Bell & de JonglbOOll) ) on a galaxy -by-galaxy basis, and systematic errors 
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non related to the choice of an universally-applied stellar IMF are 0. 1 dex. In addi tion, M/L ratios 
in sta rbursting galaxies can be biased to produce unrealistic high stellar masses (|Bell & de Jong 



20011). Those effects would have certain impact in the calculation of the SSFR, and thus, in the 



enhancement of that quantity. In order to estimate how those mass uncertainties affect our results, 
we have run two additional Monte Carlo shufflings. In Fig. [8] we show the result of this exercise. 
We have randomly added a gaussian error with a = 0.1 dex to the stellar masses of all galaxies 
and repeated the process 500 times, finding an average output value similar to the one presented 
in Fig. |4]but with larger errors. The impact of the new errors on the average enhancement (taking 
into account also the enhancement in the SFR, as we did in the previous section) is negligible. 




Fig. 8. — Enhancement in the SSFR including estimates for the errors in the stellar masses. Dotted 
line: Gaussian error with a = 0. 1 dex. Dash-dotted line: Non star forming galaxies with gaussian 
error with a = 0. 1 dex and starburst galaxies ( SFR > 16MQyr~^) with errors following an inverted 
lognormal distribution to produce a tail to the lower masses, with a shift of 0. 1 dex also to the 
lower masses and a = 0.2 dex. Solid line: Enhancement in the SSFR as in Fig!^ for comparison. 



In order to estimate the uncertainties introduced by systematics in the M/L ratio of starbursts, 
we have performed a similar exercise but using an error which includes a systematic shift down of 
0.1 dex, cr = 0.2 dex and a tail to the lower masses defined by an inverted lognormal distribution. 
We have applied this new error to objects with SFR > \6MQyr~^, which roughly corresponds to 
twice the average SFR in our sample of star forming galaxies, and the symmetric error described 
above to galaxies with SFR < l6MQyr~K The result (dash-dotted in Fig. [8]) shows some extra 
enhancement in this case, which leads to an average enhancement in the SF activity e = 1.85 ±0.35, 
barely changing the result already found. 
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Another potential source of uncertainty is the stellar masses of pairs of galaxies not resolved 
in the ground-based photometry catalog (i.e., our very close pairs group). In order to test the impact 
of this underdeblending on the galaxy masses, we take U and V rest-frame fluxes of galaxies widely 
separated, add the m together and check what stellar mass would result in the case of applying the 
Borch et aP (|2006 ) method to a galaxy with exactly the same color as the combination of the two 



galaxies, and compare with the sum of the two original masses. We find that for pairs of galaxies 
of all kinds (all-all, SF-SF and SF-all) there is a < 0.01 dex offset and 0.08 dex scatter between 
the two sets of masses. I.e., masses from combined luminosities are the same as the sum of the 
individual masses to 0.08 dex, what means that our stellar masses are extremely robust against 
underdeblending issues. 

Together with the stellar-masses, the main source of uncertainty in our an alysis is the conver- 



sion between observed 24fim and TIR in the process of obtaining the SFRs. IZheng et al.l (|2007|) 
have demonstrated that the Sbc template used here is an appropriate choice for this dataset at all 
IR luminosities, but in order to find an absolute upper limit for the final results we will show in 
^5.21 we estimate the different results we would obtain if considering an Arp220 template in some 
cases. 




Fig. 9. — Impact on the SSFR enhancement when using an Arp220 template in the conversion be- 
tween observed 24 fim and TIR luminosity for certain objects. Red line: Extreme case in which we 
apply an Arp220 template to all interacting systems (and Sbc template to everything else). Green 
line: Arp220 template applied to objects with SFR > l6MQyr~^ (and Sbc template to everything 
else). Black line: Same as in Fig. |?] for comparison. Red and green lines include the asymmetric 
stellar mass errors applied in Fig. \8i 
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We find that at a given 24/im flux, the use of an Arp220 template gives a TIR luminosity 
which is higher than that derived using a Sbc template by a factor of 2. We apply this factor of 
2 correction to the TIR luminosity of all the galaxies that we define as starburst for this purpose 
(SFR > 16M0jr"') and show the result as the green line in Fig. [9l The enhancement found in this 
case is e = 2.1 ±0.4, consistent with, by higher than the e = 1.8 ±0.3 found in ^4.1.1[ 

We also want to test the extreme case in which the IR SED of all galaxies undergoing an 
interaction follows an Arp220 SED, independently of their level of SFR. This is clearly an unre- 
alistic case as we know that some of our galaxies in close pairs and remnants have SFRs as low 
as 4-5MQyr~^ (factor of 10 less SFR than Arp220), and we know also that the averag e IR SED of 



z ^ 0.6 galaxies with SFRs > \OMQyr is similar to the Sbc template adopted here (IZheng et al 



2007bt) . but it is useful in the sense that it provides a strong upper limit beyond the uncertainties of 
our data and method. The result found in this case can be seen as the red line in Fig. |9l Clearly, 
a much stronger enhancement is present as a consequence of this overestimation of the TIR lumi- 
nosity, that leads, when taken together with the SFR enhancement calculated in the same way, to 
e = 3.1±0.6. 



4.2. How important are mergers in triggering dust-obscured starbursts? 

We have demonstrated that when averaged over all events and all event phases there is a 
relatively modest SFR enhancement from major galaxy merging and interactions. It is of interest to 
constrain how the distribution of SFRs differs between the non-interacting and interacting galaxies. 
Here we present a preliminary result on one aspect of the issue, namely the fraction of infrared- 
luminous galaxies that are in close pairs < 40kpc or were visually classified as merging systems. 

In Fig. [TO]we show how the fraction of galaxies that are either in close pairs (r^ < 40 kpc) or in 
morphologically-classified merger remnants varies as a function of their total IR luminosity. This 
fraction is constant (~20%) for 6 x lO^^L© < Ltir < 3 x IO^'Lq. At higher luminosities, the merger 
fraction increases as a function of the IR luminosity, reaching 55% just below Lrm = IO^^Lq. The 
lower IR limit of 6 x was chosen to ensure a flux of 83/iJy over the entire redshift range. 

The increase in merger fraction at high IR luminosit y is in accord with previous results at both low 



and intermediate redshift (e.g.. ISanders et al.lll988h . This suggests that merging and interactions 



are an important trigger of intense, dust-obscured star formation. Apparently, high IR luminosities 
are difficult to reach without an interaction. 

A key point, however, is that not all mergers have high IR luminosity. While mergers can 
produce enormous SFRs, and also Ltir > lO^^L© is best reached by merging, the typical SFR 
enhancement in mergers is modest. 



10" 

TIR (U) 



Fig. 10. — The fraction of systems in close projected pairs rp < 40 kpc or in visually identified 
mergers as function of the total IR luminosity of all 24iim detected galaxies in the sample. The 
merger fraction is ~ 20% from our luminosity limit of 6 x lO'^L© to 3 x IO^^Lq. Higher than this 
luminosity, the merger fraction begins to grow to 55% just below lO'^L©. 



5. Discussion 



We have assembled a unique data set for galaxies at 0.4 < z < 0.8 that combines redshifts, 
stellar masses, SFRs and HST morphologies to explore the role of major mergers and interactions 
in boosting the SFR. In practice, we have combined projected correlation-function and morpho- 
logical techniques to estimate the average enhancement of star formation in star-forming galaxies 
with M^, > 1O'"M0 and 0.4 < z < 0.8, where the average is taken over most merging phases and 
all mergers. We find a SF enhancement by a modest factor of ~ 1.8 for separations of < 40 kpc 
in both the SFR and the SSFR. How does this compare with previous observations and models? 
What implications does this mild enhancement have on the contribution of major mergers to the 
cosmic SF history? 



5.1. Comparison with previous observations 



Our analysis is most directly c omparable to estimates of star formation rate enhancement in 
galaxy close pairs by lLi et al.l(|2008|) using the SDSS at z ~ 0. 1 because of the similarities between 
our methods. Using a cross-correlation between star-forming and all galaxies, they found an en- 
hancement of ~ 1.45 for an average galaxy mass of {log(M^/MQ)) = 10.6 within a radius of 15 to 
~ 35 -40 kpc. In Fig. [TTl we show the comparison between their present-epoch measurements and 
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Fig. 1 1 . — Comparison be tween the enhancement found in this work (black filled points) and the 
one found in [Li et al\ d200q) a? z ~ 0. 1 (open diamonds). In both cases, a cross-correlation SF-all 
is shown. Both works show a statistically significant enhancement of the SSFR at rp < AQkpc. 



ours (average galaxy mass ~ IO^^^Mq, SF-all cross -correlation) revealing reasonable quantitative 
agreement. Both the projected separation scale (< 40kpc) and the overall amplitude at small pro- 
jections (x 1.5-2) agree. The enhance ment found here a lso agrees (given the error bars) with the 
enhancement found at 0.75 < z < 1.1 in lLin et al.l (l2007h . Our results are similar to those at both 
z = 0. 1 and at z = 1 , despite the factors of sev eral difference betwe en the typical star formation rates 
of galaxies between z = 1 and z = (see, e.g. JZheng et al.ll2007b|) . This is interesting, and points to 
a picture in which at least the average enhancement of star formation in galaxy interactions appears 
to be independent of the 'pre-existing' star formation in the population. 



Lin et al.l (|2007h also measured an enhancement in the TIR emission in galaxy pairs and merg- 



ers in the 0.4 < z < 0.75 range. They find that the infrared luminosity of close pairs with both mem- 
bers selected to be blue is 1.8 ±0.4 times that of control pairs, similar to our value 1.75 ±0.18 for 
close pairs in the bin ISkpc < rp < AOkpc from the SF-SF autocorrelation. For late-phase mergers, 
they measure 2. 1 ± 0.4, marginally consistent with our 1 .54 ± 0.08 from the SF-SF autocorrelation 
at rp = Okpc. Slight differences in that number may be attributed to differences in the 'merger' 
classification. For example, many of the remnants that we include in this study may not be de- 
tected by automated meth ods based on the intensity- weighted Gini-M2n or asymmetry parameters. 
As shown in recent work Jjogee et al.ll2009l : iMiller et al.ll2008|) . automated methods based on CAS 
asymmetry parameters tend to capture only a fraction (typically 50% to 70%) of the visually- 
identified merger remnants and often pick up a dominant number of non-interacting galaxies that 
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have small-scale asymmetries associated with dust and star formation. 



In related work ( Jogee et al. 2009 1). we recently estimated the overall merging rate and also 
addressed the SFR enhancement at 0.24 < z < 0.8. For the subsample of systems with > 
2.5 X lO'°A/0 we find that the average SFR of late stage mergers with mass ratio between 1:1 and 
1:10 (both major and minor mergers) are only enhanced by a modest factor (1.5-2 from their Fig. 
15) with respect to non-interacting galaxies. There are three differences that make it difficult to 
perform an exact comparison between our works: a) In the pres ent study we try to isolate the 
contribution from major interactions (mass ratio 1:1 to 1:4) while in I Jogee et al. J2OO9I) we focused 
on both major and minor interactions (mass ratio from 1 : 1 to 1 : 10). b) The normalization is slightly 
different because in t he present w o rk we compare the SFR in mergers with the SFR in the pair of 
progenitors, while in Ijogee et all toO^ we compare with the SFR of individual galaxies with 
mass similar to the interacting system. As the average SFR is a function of the galaxy stellar 



\0^^ .progenitor 



^SFR, 



M>2x 10'" descendant ■ 



c) In the present pap er we attempt to target 



mass, 2xSFRm> 

both early and late phases of the interaction, while in iJogee et al.l (|2009|) we focus on the later 
phase. Nonetheless, it is encouraging that the two studies agree qualitatively in finding a modest 
enhancement in the average SFR in galaxy interactions. 

Taken together, we argue that our results are consistent with those of previous works. We have 
used a bigger sample of galaxies with both HST/ACS and Spitzer/MIFS coverage than previous 
works at z > 0.4 and we tried to trace the SFR enhancement in all the stages of the interaction 
with a consistent treatment of ground-based selected galaxy pairs and morphologically-selected 
pairs and remnants. We view it as extremely encouraging that where the wor ks are the mos t ro- 
bust (close pa irs), the results are highly consistent (comparing our work with iLi et al.ll2008l and 
the pairs from lLin et al.l 120071) . It is clear that robustly assessing the star formation enhancement 
in advanced- stage mergers, identifiable using only high-resolution data and morphological tech- 
niques, is considerably more challenging. The results for advanced- stage mergers are therefore 
less well-constrained, but are nonetheless all consistent with a modest but significant enhancement 
in SFR. 



5.2. What fraction of star formation is triggered by major interactions? 

We can now combine our estimates for the SFR enhancement, the fraction of galaxies in 
projected close pairs, the average SFR, and the amount of SFR in recognizable merger remnants to 
quantify what fraction of star formation at 0.4 < z < 0.8 is directly triggeredhy major interactions. 
We will not include systematics such as the uncertainty in conversion of 24^m to total IR, or the 
effect of the 24^m flux limit, but we will consider the systematics driving the difference between 
the SSFR enhancement and the SFR enhancement. We make the cross-correlation between star 
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forming galaxies as primaries and all galaxies as secondaries our default case because it includes 
the residual SFR in red galaxies and also traces the SF enhancement in disk galaxies during the 
encounter with a non star forming galaxy. We will also show the values obtained for the SF-SF 
autocorrelation. We use the values 1.80ib0.30 and 1.50±0.25 found in § |4.1.1| in rp < 40kpc 
systems, for the cross -correlation and the autocorrelation respectively. 

The fraction of galax ies in close physical pairs within a separation r /■ can be d erived using the 
following approximation JPatton et allboool : iMasjedi et allbooel : IbcU et ahlbood) : 



P{r < Vf) : 



Atth 
3^ 



(7) 



Here P{r < Vf) is the fraction of galaxies in the parent sample in pairs with real separations of 
r < Tf, n is the number density of galaxies satisfying the pair selection criteria, and ro and 7 are 
the parameters of the power-law real-space correlation function of the parent sample, subjected to 
the pair selection criteria (i.e., we use a stellar mass ratio of 1:1 to 1:4 as a requirement for a pair 
to enter into the correlation function). 

Note that because in this paper we typically impose criteria for matching and forming pairs 
(e.g., a mass ratio between 1:1 and 4:1), the number density n used is not the number density of 
the larger parent sample ^parent- The number of possible pairs at any projected separation range 
is lower than in the case in which no mass ratio criteria is imposed because many pairs with 
mass ratios beyond the allowed limit are automatically rejected. As the fraction of galaxies in 
close physical pairs is directly related to the number density of galaxies n, this parameter has 
to be fine-tuned in order to get the right fraction. The number density used in Eq. |7] has to be 
corrected for the effect that the mass ratio criteria introduces on the total number of potential 
pairs. Then, the number density of the larger parent sample n^saenx is not used here, instead we use 
n = «parentA'pairs/[0.5A/parent(A'parent " 1)], whcrc A'pairs Is the numbcr of pairs that can be formed in the 
parent sample given the matching criteria, and Apparent is the number of galaxies in the parent sample 
{N{N- l)/2 is the expression for the number of possible pairs in the case of simply pairing up the 
parent sample). 



We tested this approximation using the semi-analytic galaxy catalog of lDe Lucia et al.l(|2006h . 
derived from the Millennium A^-body simulation. At z ~ 0.6 these simulations matched well the 
stellar mass function and correlation function of > IO^^Mq galaxies. We find that at r < 50kpc 
Equation |7] is a good approximation to the actual fraction of galaxies in close pairs in the simu- 
lation; at larger separation Equation [7] is increasingly incorrect (this is the subject of a paper in 
preparation). 



From fits to the projected two-point cross-correlation function of our sample, we determine 
ro = 1.8 ± 0.2 Mpc, 7 = 2.2 ± 0. 1 for the real-space correlation function, and n = 0.0152 galaxies 
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per cubic the latter gives a sample of A'gai = nV = 1913 galaxies in the volume probed by this 

study. This yields P{r < 40kpc) = 0.06 ±0.01 (i.e., 6% of sample galaxies are in close pairs with 
real-space separations < 40kpc). With this real-space two-point correlatio n function, 75± 10% 



of all projected close pairs should be real close physical pairo (Eq. 6 of IBell et al.l 120061 and 
confirmed using the Millennium Simulation at the redshift of interest). Thus the fraction of objects 
in projected close pairs will be /pair.proj = 0.06/0.75 = 0.08 ±0.02. This fraction includes projections 
due to real structures like clusters or filaments, but not the purely random projections due to redshift 
uncertainties which would be present if we would just count the galaxies in projected close pairs 
in our catalogue. 

When considering all pairs at all separations in our sample with > lO'^M©, mass ratios 
between 1 : 1 and 1 : 4, and primary galaxies with lAfim fluxes > 83/i Jy and/or blue, the average 
SFRis (5Fi?)jypi^^ip^i^ = 13.2±0.6MQyr-^ The toto/SFR in the M.em = 38 ±5 recognizable merger 
remnants is 5Fi?remnants = 753 ±97A/0yr"^ Thus, we can calculate the fraction of SFR occuring in 
pairs with separations < 40kpc: 

■^gal/pair,proj0.5e(5Fi?) .^j . -|-5Fi?i.emnants 

:; — : = ZU±jyo (Oj 



A^gaiO.5 {SFR) 



typical, pair 



for the SF-AU correlation. A similar analysis with the SF-SF correlation yields 16 ±3%. What 
we have done in the numerator of Eq.[8]is to take the typical SFR in our pairs and divide by two in 
order to get the typical SFR of a galaxy contributing to such pairs. This number is different from 
the typical SFR in our galaxy sample for two reasons: first, we have imposed a mass-ratio criterion 
(only allow pairs with mass ratios between 1 : 1 and 1 :4), which makes the averaged SFR in all pairs 
to be slightly biased high respect random pairs without any mass ratio criterion, and second, the 
fact that we force the primary galaxy to be a star former (in the case of the cross -correlation) has 
a similar effect. Then we have multiplied it by the enhancement e in order to take into account 
the excess SFR triggered by major interactions and introduced the factor Ngai to account for all the 
SFR occuring in those galaxies. A key piece of Eq.[8]is the different treatment of merger remnants. 
The correlation function can tell us what is the fraction with separations between rp = 40 kpc and 
rp = kpc but we have defined merger remnants as objects which have already coalesced, so if we 
think in terms of the duration of the interaction instead of the separation between the galaxies, these 
objects would be beyond the reach of the correlation function, and have to be treated separately. In 
the denominator we have only divided by the total SFR occuring in all the galaxies contributing 



^The correlation function is calculated in proper coordinates, because the process of interest is galaxy merging and 
close pairs of galaxies have completely decoupled from the Hubble flow. 

This is only valid after removing the effect introduced by purely random projections with the correlation function 
method. 



-32- 



to any pair we can fomi with the already mentioned criteria. The difference between this factor 
and the total SFR calculated simply adding up the SFR of all galaxies in the sample is 5% and is a 
consequence of the few galaxies which are not paired with any other galaxy in the sample. 

Yet, the fraction of the total SFR that occurs in pairs and remnants with < 40kpc separation 
does not immediately characterize the SFR triggered by interactions, because ~ 12% of SF should 
happen at rp < 40kpc anyway, as we show below. Only the excess star formation in pairs and 
remnants should be atributed to triggering by interactions: 



A^gal/pair,projO-5(e- 1) (5Fi?)typical,pair + ('^^^remnants " A^mn ('^^^)typical,pair) _ o _^ 3<v 



(9) 



Again, a similar analysis for the SF-SF autocorrelation yields 5 ± 3%. These values for the excess 
are 12% lower than those in Eq.[8]due to the total number of interacting systems, which is higher 
than the 8% of galaxies in close pairs mentioned before because it includes the merger remnants 
that are not taken into account by the correlation function method. 

Taking all this together, this analysis shows that only ~ 8% of the star formation at 0.4 < 
z < 0.8 is triggered by major mergers/interactions. This may seem in disagreement wit h previous 



results fro m 'morphological' studies. We therefore compare our results with those of iBell et al 
hoosi and lWolf etaP fcoOsL who found that ~ 30% of the globa l SFR at z = 0.7 is taking place 



in morphologically perturbed systems and with ljogee et al.l (|2009h . where we find a similar result 
at0.24<z<0.8. 



Both lBelletal.1 (l2005h and lWolf et all (l2005h performed a study of the total SFR occuring in 
visually-classified interacting galaxies in a thin redshift slice 0.65 < z < 0.75 without imposing 
a lower mass limit (only an apparent magnitude limit). That is the key difference between those 
earlier works and ours. We impose a mass cut in this paper of IO^^Mq and 2 x IO'^Mq for galaxies 
and visually-classifi ed interact i ons re spect ively. E.g. t he fra ction of SFR in galaxies with > 
1 - 2 X lO^^M© that IBell et all J2OO5I) and Iwolf et aP J2OO5I) identified as interacting/peculiar is 
15-21%, co mpared to 20% in Eg. [81 Only 1/6 of the star formation in the interacting/peculiar 
galaxies from lBell et al.l (|2005|) and lWolf et al.l (|2005[) occurs in what we would designate as merger 
remnants, with the other 5/6 occuring in galaxy pairs. Jogee et al. (2009) argued that 30% of star 
formation was in systems that they classified as major or minor interactions, with mass limits 
different from those used in this paper. This value is an upper limit to the fraction of star formation 
in major mergers where each galaxy has mass > IO^^Mq, both because of the effect of mass limits, 
and because minor mergers host much of the star formation in systems that they classified to be 
interacting. Accounting for these differences, our result is in qualitative agreement with theirs. 



However, the key difference is that neither lBell et al.l(l2005|) . IWolf et al.l(l2005|) nor ljogee et al 



(|2009|) try to quantify the excess of SF in interacting systems, as we do in going from Eq. [8] to 
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Eq.[9l In summary, our new results here pose no inconsistency with earlier studies, but refine them 
by quantifying the physically more relevant quantity of SF excess. 

We have presented the average enhancement in SFR caused by major mergers of galaxies with 
masses above IO^^Mq at 0.4 < z < 0.8, deriving that approximately 8% of the SF in the volume is 
directly triggered by major merging. As we have mentioned before, the SFR enhancement e seems 
to be roughly independent of the quiescent SFR ground level present in the galaxy population, 
that is, insensitive to the drop in the SFR density of the Universe since z=l. If this is true, it 
means that the fraction of star formation directly triggered by galaxy interactions (given a mass 
cut) would depend only on the number o f galax ies undergoing interactions. Using the evolution 



of pair fraction found in iKartaltepe et al.l (|2007|) we can infer a directly-triggered SF fraction of 



1-2% in the local Universe, as well as a fraction of 14-18% at z=l. On the other hand, assuming 
no evolution in the pair fraction would keep the merger-triggered fraction at 8% between z = and 
z = 1 . These numbers have to be taken extremely carefully by the reader, as we present here only a 
crude extrapolation of our results to different redshifts in order to get an idea of the importance of 
the merger-driven star formation in the Universe. 

There are a number of limitations of our result that should be borne in mind. First, we can 
only include star formation rates > 5Moyr"^ for z ~ 0.6 galaxies. Therefore, our estimates of the 
SF contribution from merging may be an upper limit because the merger-driven boost in SFR will 
cause more objects to satisfy this criterion. Second, there are uncertainties in the conversion of 
24/ir n to total IR, which coiild influence the ex cess star formation in close pairs or mergers by ~ 



30% (|Papovich & Bellll2002l : IZheng et al.ll2007|) : this could be addressed once longer-wavelength, 
deep Herschel RA.CS observations become available. We can only calculate an absolute upper 
limit by using the value of the enhancement found in the extreme case in which the 24/im to TIR 
conversion in all the interactions, independently of their luminosity, is calculated using an Arp220 
template. Using as input e = 3. 1 ± 0.6 for Eq. |9] we would find a directly triggered fraction of 19 ± 
5%, consistent with an scenario in which the underlying level of SF is basically negligible and most 
of the new stars are being formed in the burst mode. Third, it is conceivable that some enhanced star 
formation occurs in very late- stage merger remnants that were no longer recognized as remnants 
and hence were not included in this census. This is both a practical (classification) and conceptual 
issue: when does one declare a merger remnant a normal galaxy again? Nonetheless, despite these 
points, the analysis presented here has made it very clear that only a small fraction of star formation 
in galaxies with M* > 10 "'M© at0.4<z<0.8is triggered by major interactions/mergers. 



5.3. Comparison with theoretical expectations 



Star formation enh ancement in mergers has been studied extensively with hydrodynamical N- 
body simulations (e.g Isarnes & Her nQuist*^199ll: iMihos & Hemquistl[l994l 1 19961 : ISpringellllOOol : 
Cox et al. Il2006a[ l2008hldi Matteo et al. 200 7). However, large-scale cosmological simulations lack 
the dynamic range to resolve the internal dynamics of galaxies, crucial for modeling the gas inflows 
and the associated enhancement in star formation. Therefore, the majority of these studies (except 
Tissera et al.ll2002l) have been of binary galaxy mergers with idealized initial conditions, typically 
bulgeless or late-type disks. In most studies, the properties of these progenitor disks are chosen to 
be representative of present-day, relatively massive spiral galaxies such as the Milky Way. These 
studies have shown that the burst efficiency in mergers is sensitive to parameters such as merger 
mass ratio and orbit, and progenitor gas fraction and bulge content. Therefore, any attempt to 
use these results in an ensemble comparison must somehow convolve these dependencies with a 
redshift dependen t, cosmologically motivated distribution function for these quantities. In addition. 
Cox et al.l (|2006ah have shown that star formation enhancement in mergers can also depend on 
the treatment of supernova feedback in the simulations. Furthermore, the detailed star formation 
history during the course of a merger, pa rticularly in the late s tages, may depend on the presence 
of an accreting supermassive black hole (Idi Matteo et al.ll2005|) . 



Let us conside r the results from representative examples of such binary merger simulations, 
by lCox et al.l (|2008|) . who studied a broad range of merger mass ratios, gas fractions, and progeni- 
tor B/T ratios, as well as exploring the effects of two different SN feedback recipes. The 1 : 1 merger 
of two "Milky Way"-like progenitors (shown in their Figure 12) shows an average factor of ~ 1.5 
enhancement in SF over about 2.5 Gyr, and a larger enhancement of a factor of 2-10 for a shorter 
period of about 0.6 Gyr. The overall average enhancement over the whole merger is about a factor 
of 2.5, depending on the precise timescale one averages over. Very large enhancements (~ 5-10) 
occur over a very short timescale, < 100 Myr. This particular simulation represents the largest 
expected SF enhancement, as the burst efficiency incr eases strongly to wards equal merger mass 
ratio. For mergers with 1:2.3 mass ratio (Figure 10 of ICox et al.1120081) . there is an enhancement 
of a factor of ~ 1.5 for 2.5-3 Gyr, and of 2.5 for about 0.6 Gyr. It is also interesting to note that 
the star formation rate in the l ate stages of the merger, when the galaxy still appears morpholog- 
ically disturbed (see Fig. 7 of ICox et al.ll2008|) is depressed with respect to the isolated case. A 
di verse set of progenitor morphologies, ranging from ellipticals to late-type spirals, w as studied 



by Idi Matteo et al. J2OO7I) . Overall, their results are qualitatively similar to those from ICox et alj 
J2OO8I) 



The simulations discussed so far aimed to reflect progenitor disks with gas fractions, sizes, 
and morphologies typical of relatively massive, low-redshift late-type spir als such as the Milky 
Way: gas fraction fg ~ 0.2; B/T ~ 0.2; scale length ~ 3 kpc. However. iHopkins et al. (l2009i) 
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show that the burst efficiency is strong function of progenitor gas fraction, in the sense that higher 
gas fraction progenitors have weaker fractional enhancements. The burst efficiency is a factor of 
eight lower for a gas fraction of 90 % than for the canonically-used value of 20%. It is worth noting 
here that we find the same level of SF enhancement in major mergers at z ~ 0.6 and z ~ 0. 1, where 
the gas fractions of the two samples are expected to be rather different (se e S5.ll) . Whether or not 
this is quantitatively at odds with t he expectations of [Hopkin s et al.l (|2009l) remains to be seen. On 
the other hand, recent results from ldi Matteo et al.l (1200 8) show no difference between the strength 
or duration of tidally-triggered bursts of star formation in local Universe and their higher redshifts 
counterparts, in good agreement with the present study. 



To place results in a cosmological context. ISomervil 



large suite of hydrodynamic merger simulations (|Cox et al 



eetal 



2006b 



(l2008h used the results from a 



Robertson et al.ll2006alJb| ]c) to 



parameterize the dependence of burst efficiency and timescale on merger mass ratio, gas fraction, 
progenitor circular velocity, redshift, and the assumed effective equation of state. They imple- 
mented these scalings within a cosmolo gical semi-analyti c merg er tree model. We applied our 
selection criteria to mock catalogs from ISomerville et al.l (|2008|) . by comparing the fraction of 
SFR produced in the triggered mode in galaxies with > 2 x lO'^M© in our redshift range which 
suffered a major merger in the last 500 Myr with the total SFR occuring in galaxies > IO'^Mq. 
We found that approximately 7% of the SFR in the volume is produced in the burst mode triggered 
by major mergers. This is in excellent agreement with the 8 ± 3% of the overall SFR being directly 
triggered by major interactions we showed in the previous section. 



6. Conclusions 

To quantify the average effect of major mergers on SFRs in galaxies, we have studied the 
enhancement of SF caused by major mergers between galaxies with > IO^^Mq at 0.4 < z < 0.8. 
We combined redshifts and stellar masses from COMBO-17 with high-resolution imaging based 
on HST/ACS data for two fields (ECDFS/GEMS and A901/STAGES) and with star formation 
rates that draw on UV and deep 24/zm data from Spitzer to form a sample a factor of two larger 
than previous studies in this redshift range. We then applied robust two-point correlation function 
techniques, supplemented by morphologically-classified very close pairs and merger remnants to 
identify interacting galaxies. Our main findings are as follows: 

1 . Major mergers and interactions between star- forming massive galaxies trigger, on average, a 
mild enhancement in the SFR in pairs separated by projected distances rp < 40kpc; we find 
an enhancement of e = 1 .80 ± 0.30 considering the SF-AU cross-correlation, where only one 
galaxy in the pair is required to be forming stars. For a similar analysis using the autocorre- 
lation of star forming galaxies we find e = 1 .50 ± 0.25. 
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2. Our results agree well with previous studies of SF enhancement using close pairs at z < 1 . 
In particular, the behavior of SF enhancement at z = 0. 1 , z = 0.6 and z = 1 appear to be rather 
similar, indicating that the average SFR enhancement in galaxy interactions is independent 
of the 'pre-existing' SFR in the population. 

3. We combine our estimate of the average SFR enhancements in major mergers with the global 
SFR to show that overall, 8 ± 3% of the total star formation at these epochs is directly trig- 
gered by major interactions. We conclude that major mergers are an insignificant factor in 
stellar mass growth at z < 1 . 

4. Major interactions do, however, play a key role in triggering the most intense dust-obscured 
starbursts: We find that the majority of galaxies with IR- luminosities in excess of 3 x 10^ ^L© 
are visually classified as ongoing mergers or found in projected pairs within < 40kpc sep- 
aration. This is not in disagreement with the small average SFR enhancement if the most 
intense SF bursts last only ~ lOOMyr. 

5. Our results for the SF enhancement app ear to be in qual i tative agreem ent with the extensiv e 
suite of hydrodynamical simulations by Idi Matteo et all 120071 l2008b and lCox et all pOOSb . 
who produce both intense, short-lived bursts of SF in some interactions, but yet produce 
average enhancements of only 25-50% averaged over the ~ 2 Gyr timescale taken to com- 
plete the merger. Furthermore, we find excellent agreement between the fraction of the total 
SFR directly triggered b y major merging measur ed here and the 7% calculated from mock 
catalogues obtained from lSomerville et al.l (|2008|) . 
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